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Transverse Spectra of Hadrons in Central AA Collisions at RHIC and 
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q ■ We study the transverse spectra of hadrons in nearly central AA collisions at RHIC 
and LHC in a broad transverse momentum range pQ. Low-pr spectra are calculated by 
using boost-invariant hydrodynamics with initial energy and net-baryon densities from 
the EKRT [2] pQCD+saturation model. High-p T spectra are obtained from pQCD jet 
calculation |3] including the energy loss of the parton jlj in the matter prior to its frag- 
mentation to final hadrons. 
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Transverse momentum spectra of hadrons in ultrarelativistic nuclear collisions provide 
valuable information of particle production mechanism in the collisions as well as dynamics 
and properties of the produced QCD matter. Low-pr features of the single particle spectra 
are well described by hydrodynamical models and the data is consistent with an ideal 
fluid behaviour of the matter. The observed azimuthal asymmetry in non-central Au+Au 
collisions at RHIC has been argued to result from strong collective motion and early 
thermalization of produced partonic matter. The suppression observed in the high-p^ tail 
of the spectra relative to the p+p and d+Au collisions is understood as the energy loss 
of high-pr partons into the thermalized partonic matter. 

Our approach is to calculate initial particle production using the EKRT model which is 
based on the idea that the low-p^ particle production is controlled by saturation among 
the final state gluons in contrast to the initial state saturation models, where saturation is 
the property of colliding nuclei. The EKRT model provides a closed framework to calcu- 
late initial transverse energy and net-baryon number in midrapidity in nuclear collisions 
with sufficiently large \fs and A. Final low-pp hadron spectra are calculated by using 
hydrodynamics with initial densities from the EKRT model. A good agreement with the 
measured data in central Au+Au collisions at RHIC is obtained jlloj . We further predict 
the hadron spectra in central Pb+Pb collisions at the LHC PQ. 

High-pr part of the spectra is calculated using factorized pQCD parton model for high- 
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Pt parton production and taking into account the parton energy loss before fragmentation. 

We compare these two models with the RHIC data and determine the regions where 
each component is dominant. In particular, we find that the low-j>r hydrodynamical 
spectrum dominates over the fragmentation spectrum in a much wider px region at the 
LHC than at RHIC. We discuss the independence of the two components and conclude 
that they are most likely almost independent even in the cross-over region. 

2. Models 

2.1. pQCD + saturation + hydrodynamics 

The EKRT model [2. estimates the final state saturation scale using the following 
geometric criterion: Saturation becomes important when the produced gluons fill the 
whole available transverse area of the colliding nuclei. For central collisions this can be 
written as 



where Naa is the number of gluons above a transverse momentum cut-off, qr > qo in 
the rapidity interval \y\ < 0.5, and Ra is the nuclear radius. This condition provides the 
saturation scale p sat = qo- If p sa t 3> ^qcd, pQCD can be used to estimate the number 
of produced partons and the amount of transverse energy at midrapidity. This approach 
gives also the net-baryon number at midrapidity. If we assume that the produced matter 
thermalizes immediately after production at r = l/p sa t, we obtain the initial energy and 
net-baryon density at r for hydrodynamical evolution. 

We use ideal fluid hydrodynamics with boost-invariance and azimuthal symmetry. In 
the bag-model equation of state the hadron gas phase consists of all hadronic states with 
m < 2 GeV and the QGP phase of massless gluons and three flavors of quarks. Critical 
temperature is chosen to be 165 MeV. After the hydrodynamic expansion and cooling 
of the matter, the Cooper-Frye decoupling prescription is applied for the calculation of 
the low-pr spectra. Below we show the sensitivity of the decoupling procedure on the 
decoupling condition. 

2.2. pQCD + fragmentation + energy loss 

High-pT spectra are calculated using nuclear parton distributions, pQCD parton cross 
sections, fragmentation functions and quenching weights for energy losses. We use leading 
order perturbative cross sections with K-factors fixed from p+p(p) data The K-factors 
are extrapolated to the LHC energy by using different parametrizations to estimate the 
uncertainties in the extrapolation. The magnitude of the energy losses can be expressed 
with one effective transport coefficient, which is fixed from the RHIC Au+Au data at 
y/SNN = 200 GeV jlj. Since the transport coefficient is proportional to energy density, 
fixing it in one collision predicts it for other collisions. There is quite a large uncertainty 
associated with the eikonal approximation used in the energy loss calculation; arbitrarily 
large energy losses are allowed whereas the energy of real jets is limited. There are different 
ways to deal with this j3], shown here as an uncertainty in our pQCD fragmentation + 
energy loss results. 
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3. Results 

In Fig^i we show our results for unidentified charged hadrons at midrapidity for 5 % 
most central Au+Au collisions at ^/s^n = 200 GeV. As mentioned above we show our 
hydro results with two different decoupling temperatures to show the sensitivity on the 
decoupling condition. From the figure we observe that both the normalization and slope 
of the spectra are well reproduced with a single decoupling temperature Td ec = 150 MeV. 
The spectra from fragmentation calculation with and without energy loss are plotted 
in the same figure. It is seen clearly that the measured spectra at high px cannot be 
explained without energy loss. The transport coefficient that determines the energy loss 
is fixed from this data so that in this case the agreement is obtained by construction. 
The band in the energy loss results shows the uncertainty of eikonal approximation. For 
identified particles the spectra can be found from pQ. 

Fig shows our prediction for the 5 % most central Pb+Pb collisions at the LHC. 
The hydrodynamical results are again shown as a band corresponding to the decoupling 
temperatures between 120 and 150 MeV. The band for the pQCD fragmentation results 
without energy loss is from the uncertainty in the extrapolation of the K-factors to the 
LHC energy. The large uncertainty in the energy loss calculation comes, as before, from 
the eikonal approximation. 

We see in Fig that at RHIC the hydrodynamical and pQCD + energy loss spectra 
cross at p-r ~ 3 . . . 4 GeV, and that the data starts to deviate from the hydro spectrum 
in the same pt region. At the LHC, the crossing region is moved to px ~ 5 ... 6 GeV, 
suggesting a wider px region of applicability of hydrodynamics at the LHC than at RHIC. 
It is also interesting to note that in the cross-over region the fragmentation and the 
hydrodynamical components are most likely almost independent: At RHIC 95 % of the 
thermalized matter comes from mini-jets with partonic transverse momenta p sat < qr < 
3.6 GeV and the higher-g^ partons contribute to the normalizations and the slopes of 
the hydrodynamic spectra only slightly. On the other hand even without the energy 
loss the pQCD pions are dominantly from partons with qx ~ \HpT-, which is 5.1 GeV 
for px ~ 3 GeV pions. With energy loss included they originate from even higher qr- 
Thus the partonic origin of fragmentation pions and pions from thermalized matter is 
quite different suggesting that the two contributions are quite independent even in the 
cross-over region and can be added without serious double counting. The same argument 
holds also at the LHC. Crossing between hydro and fragmentation spectra depends also 
on particle species as studied at RHIC energies in ref. [§]. 

4. Conclusions 

We have calculated Iow-pt spectra of hadrons for central Au+Au collisions at RHIC 
and Pb+Pb collisions at the LHC using the EKRT model to calculate the initial parton 
production and hydrodynamics to calculate the expansion of produced matter. High-j>r 
hadron spectra are calculated by assuming that the high-g-r partons do not thermalize but 
fragment to hadrons after loosing energy when traversing the thermal matter. Our model 
is shown to be in a good agreement with the measured data in central collisions at RHIC. 
We have provided predictions for central Pb+Pb collisions at the LHC, including also 
the net baryon number, see pQ. The origins of the thermal and the pQCD fragmentation 
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spectra are discussed and it is argued that even in the cross-over region where the two 
components are comparable, they are essentially independent. 
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Figure 1. a) Transverse momentum spectra of charged hadrons at y = in 5 % most central 
Au+Au collisions at ^/sN~N = 200 GeV. The solid and dotted lines shows our hydrodynamic re- 
sults with Tdec = 150 MeV and Td ec = 120 MeV respectively. The shaded band shows the pQCD 
fragmentation+energy loss spectrum and dashed line pQCD fragmentation without energy loss. 
The data is from Refs. |7I8I9I10| . b) As Fig. but for the 5 % most central Pb+Pb collision 
at Js^N = 5500 GeV. 
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